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A combined experimental and theoretical investigation of the electronic structure of the archetypal 
oxide heterointerface system LaAlOa on SrTiOs is presented. High-resolution, hard x-ray photoe- 
mission is used to uncover the occupation of Ti 3d states and the relative energetic alignment - and 
hence internal electric fields - within the LaAlOs layer. Firstly, the Ti 2p core level spectra clearly 
show occupation of Ti 3d states already for two unit cells of LaAlOs. Secondly, the LaAlOs core 
levels were seen to shift to lower binding energy as the LaA103 overlayer thickness, n, was increased 
- agreeing with the expectations from the canonical electron transfer model for the emergence of 
conductivity at the interface. However, not only is the energy offset of only ~300meV between n = 2 
(insulating interface) and n — 6 (metallic interface) an order of magnitude smaller than the simple 
expectation, but it is also clearly not the sum of a series of unit-cell by unit-cell shifts within the 
LaAlC>3 block. Both of these facts argue against the simple charge-transfer picture involving a cu- 
mulative shift of the LaAlC>3 valence bands above the SrTiC>3 conduction bands, resulting in charge 
transfer only for n > 4. We discuss effects which could frustrate this elegant and simple charge 
transfer model, concluding that although it cannot be ruled out, photodoping by the x-ray beam 
is unlikely to be the cause of the observed behavior. Turning to the theoretical data, our density 
functional simulations show that the presence of oxygen vacancies at the LaA103 surface at the 25% 
level reverses the direction of the internal field in the LaAlOs. Therefore, taking the experimental 
and theoretical results together, a consistent picture emerges for real-life samples in which nature 
does not wait until n = 4 and already for n = 2, mechanisms other than internal-electric- field-driven 
electron transfer from idealized LaAlOs to near-interfacial states in the SrTiC>3 substrate are active 
in heading off the incipient polarization catastrophe that drives the physics in these systems. 



PACS numbers: 73.20.-r, 73.40.-c, 73.61.-r 
I. INTRODUCTION 



Since they were first studied in 2004 by Ohtomo and 
Hwang, SrTi0 3 /LaA10 3 (STO/LAO) heterostructures 
have attracted a great deal of interest. 1-4 The nature 
of the conduction that occurs when an LAO film four 
or more unit cells (UC) thick is grown on a Ti02 termi- 
nated STO substrate, 5 is still a subject of intense debate. 
The so-called 'polar catastrophe' is proposed to play a 
central role in the explanation of this conductivity. 1 Ac- 
cording to this model, the electric potential should in- 
crease linearly with each layer of LAO that is added 
to the sample. Because this cannot continue indefi- 
nitely, there must be a critical thickness at which the 
system finds a way of relaxing to a lower energy state. 
The most widely discussed model entails an electronic 



reconstruction 6 (electron transfer picture) whereby elec- 
trons are transferred from the outermost AIO2 layer to 
the Ti02 layer at the interface between STO and LAO. 
Theoretical calculations have predicted a shift of the 
LAO energy bands of 0.9 cV per monolayer of LAO.' 
For LAO overlayers four or more monolayers thick, the 
top of the LAO valence band should then rise above 
the bottom of the STO conduction band, leading to the 
transfer of electrons and averting a potential runaway. 
The transferred electrons (half an electron per unit cell) 
then form a conducting quasi two dimensional electron 
gas (q-2DEG). 6 Despite the simplicity and elegance of 
this model, as yet such shifts have not been experimen- 
tally observed in oxide MBE 8 or low O2 partial pres- 
sure PLD 9 grown thin films. As the diverging poten- 
tial, which would certainly exist for an ideal film stack, 
still needs to be avoided, additional forms of potential 
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compensation or reduction should be considered. One 
such possibility is the electron doping that would result 
from oxygen vacancies being formed in the overlayer 10-12 
or in the substrate. 13-16 Removing an oxygen ion from 
the crystal donates two electrons to the system. These 
electrons could then occupy the Ti 3d states which are 
the lowest energy electron-addition states of the whole 
system, moving to the STO/LAO interface and reduc- 
ing the built-in potential. Intermixing of cations around 
the interface is often mentioned as an alternative to elec- 
tronic reconstruction, 17-20 leading to the formation of 
Lai-^Sr^Ali-yTij, 20 or Sri-i^La^O interfaces. 21 Note 
that A-site intermixing by itself does not resolve the po- 
lar catastrophe, and it has been pointed out that not 
diffusive cation disordering but only external (redox) pro- 
cesses can alter the interfacial net charge in systems such 
as LAO on STO. 22 Many of these alternatives to the 
built-in potential paradigm involve energetics and kinet- 
ics that could differ from one film growth laboratory and 
technique to the next, yet the critical LAO overlayer 
thickness for metallic behaviour of n = 4 unit cells seems 
to be a remarkably robust phenomenon. Understanding 
this is one of the challenges currently facing the field and 
makes comparison between results from films grown us- 
ing different techniques very important. 

In order to fully understand these systems, knowledge 
of the conductive layer is vital. The electrostatic fields 
predicted by the electron transfer picture are so large 
that they should be easily observable as shifts and (ap- 
parent) broadening of LAO core levels on increasing the 
thickness of the LAO overlayer. 8 What is measured in 
a core level experiment is then a weighted total average 
of core level peaks from each LAO layer. As a result, 
a potential build-up of 3.2 eV for - for example - 4 lay- 
ers of LAO would result in significantly broadened La 
and Al core levels whose maxima would then be shifted 
by ~ 1.9 eV, due to the integration in the experiment 
over the different layer-by-layer shifts. Hard x-ray pho- 
toemission spectroscopy (HAXPES) is a powerful tool 
in detecting such core level shifts in buried interfaces. 23 
Here we present the results of HAXPES experiments on 
PLD-grown STO/LAO heterointerfaces. We turn our at- 
tention in particular to how the LAO core levels behave 
as the LAO overlayer thickness increases. We present 
non-zero LAO core level shifts of ~ 300 meV between 
two and six layers of LAO. The experimental data show 
a built-in potential of the same sign as predicted by sim- 
ple polarisation models. Its magnitude, however, is one 
order lower, and not of a unit-cell by unit-cell type. We 
present a DFT-based model in which 25% of the oxygen 
atoms in the top AIO2 layer have been removed. The 
reduction of potential build-up by oxygen vacancies is 
shown to be more than sufficient to account for the only 
modest electrostatic shifts observed in experiment, and 
would be in keeping with our experimental observation 
of Ti 3+ centers already for only n — 2. We also discuss 
alternative scenarios, including whether the experiment 
itself interferes with the build-up or maintenance of the 



proposed internal field, and although photodoping can- 
not be ruled out, we present data implying this scenario 
is unlikely. 



II. EXPERIMENT 

The measurements were performed on STO(OOl) sub- 
strates with LAO overlayer thicknesses 2 < n < 6 UC, 
where n is the number of LAO layers. In order to in- 
vestigate the effect of thick LAO overlayers an addi- 
tional n = 20 UC sample was grown. The samples are 
grown by pulsed laser deposition in an oxygen partial 
pressure of 2 x 10 -3 mbar. Oxygen partial pressure is 
known to strongly affect the transport properties of these 
systems. 24 At this relatively high pressure the presence 
of oxygen vacancies is believed to be strongly reduced 
compared to synthesis at lower partial pressures. Details 
of the sample growth have been published elsewhere. 2 

The experiments were performed at the double crys- 
tal monochromatized KMC-1 beamline 25 using the HiKE 
endstation 26 at Helmholtz Zentrum Berlin. Unless stated 
otherwise, a photon energy of hv = 3 keV was used, and 
the total energy resolution was set to be ~ 500 meV. 27 No 
surface treatment was necessary because of the large elec- 
tron inelastic mean free path A at this photon energy (~ 4 
nm, as calculated by the TPP-2M relation. 28 ) Additional 
experiments were carried out using a laboratory Al: K Q 
source (I486 eV, A = 1.9 nm) and with synchrotron radi- 
ation of energies spanning 2.15 keV to 6 keV (A = 2.8 nm 
to A = 7.4 nm, respectively). All experiments were car- 
ried out at room temperature and all energy referencing 
was done by measuring the energy position of Au 4f core 
levels. Contact between the sample and the gold refer- 
ence was made using standard silver paint, and - apart 
from a 2 mm bare space - the film surface was coated with 
commercial colloidal graphite contacting agent. During 
the photoemission experiments themselves - both for the 
synchrotron-based as well as Al: K Q -based - we checked 
for photon flux dependent effects by varying the pho- 
ton flux over more than an order of magnitude, confirm- 
ing that the measured kinetic energy, peak positions and 
widths remained unchanged. In order to eliminate even 
the smallest error caused by residual electrical charging 
in the La 4d and Al 2s core level measurements, we deter- 
mine the difference in binding energy between these core 
levels and Sr 3d for each sample and use this as our main 
quantity of interest. An inelastic electron background 
has been subtracted from the data, and they are normal- 
ized to have equal area over each of the binding energy 
regions shown. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

In Fig. la the Ti 2p spectra are shown for a bare STO 
substrate and the n = 2 and 6 samples. As in Ref. 29 
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FIG. 1. (color online) a) The Ti 2p spectrum for n — 2 and 
n — 6. As a reference the spectrum for pure STO (n = 0) 
is also included. The shoulder (inset) at 457 eV shows the 
existence of Ti 3+ (3d 1 entities) for samples both above and 
below the critical thickness of four layers of LAO. The n = 2 
and the n = 6 spectra have been shifted to lower binding 
energy (160 and 220 meV, respectively) in order to align the 
main line easing comparison between the shoulders, b) The 
Ti 2p spectrum for the n — 6 sample as a function of emission 
angle 6 between 5° and 50° away from normal emission (NE). 
Within the noise the shoulder (inset) shows marginal or no 
angular dependence, c) A simulation - based on the model of 
Ref. 29 - of the expected Ti 3+ shoulder intensity for several 
degrees of spatial confinement of the q-2DEG for an electron 
emission angle of 50° from normal emission (NE). 



all STO/LAO samples show a shoulder at the low bind- 
ing energy side of the Ti 2p 3 / 2 peak (in our data at 2.3 
eV lower binding energy), which is commonly accepted 
to be a signature of partial occupation of Ti 3d states 
i.e. the presence of Ti 3+ ions. 30 In order to facilitate 
comparison of the shoulders in the figure, the n = 2 and 
n = 6 spectra have been shifted to lower binding energy 
by 160, and 220 meV respectively to coincide with the 
STO spectrum. The observed absolute Ti 2p 3 / 2 bind- 
ing energy differences for the different LAO thicknesses 
were similar to those observed in Ref. 31. The values 
of the binding energies of the peak maxima of the Ti 
2p 3 / 2 and the Sr 3ds/2 core levels along with their ob- 
served full-width at half-maximum (FWHM) values are 
given as a function of n in Table I. We also include the 
binding energies and FWHM of the relevant LAO core 
levels. Without going into a blow-by-blow account of all 
levels for all values of n, we do mention that the Ti (Sr) 
core levels show a significant shift of 160 (240) meV to 
higher binding energy on addition of LAO on top of the 
STO substrate. The subsequent increase in the binding 
energy of these levels as n increases is limited to 70(-20) 
meV for the Ti 2p (Sr 3d) levels. 



In Fig. lb, we show that for an n = 6 film there is no sig- 
nificant dependence of the strength of the Ti 3+ shoulder 
on the electron emission angle for angles between 5° and 
50° away from normal emission (NE). Similar measure- 
ments on samples with other LAO thicknesses yielded the 
same result. 

We simulate the expected Ti 3+ shoulder using the 
model adopted in Ref. 29: 
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Here 8 is the emission angle of the excited photoelectrons, 
p is the fraction of Ti 3+ ions in the q-2DEG per unit cell, 
d is the thickness of the q-2DEG and A is the inelastic 
mean free path length of 4 nm for the 3 keV photons 
used. In Fig. lc we show a simulation of the Ti 3+ shoul- 
der for an n = 6 sample for four different thicknesses of 
the q-2DEG, 9 = 50° and p = 0.03. Comparison of the 
simulation with the experimental data (see inset) shows 
that the the confinement of the 3d 1 Ti ions (Ti 3+ ) is cer- 
tainly to a region greater than 10 UC, and more likely 
to 20 UC (8 nm) or more from the interface. This value 
is a factor of two larger than the upper limit reported 
in Ref. 29 for confinement lengths in the n = 6 samples 
studied there. Thus, the main message at this stage is 
that the Ti 3+ ions picked up in our core level photoemis- 
sion experiments are certainly not tightly confined in a 
2D layer, for example right at the STO/LAO interface. 

The second clear message our Ti 2p data convey is that 
the Ti 3d states are partially occupied already at n = 2 
(see Fig. la), even though conductivity in transport for 
the same samples only sets in at n — 4. The origin of 
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STO core levels 

n Ti 2p 3/2 BE Ti 2p 3/2 FWHM Sr 3d 5/2 BE Sr 3d 5/2 FWHM 





(eV) 


(eV) 


(eV) 


(eV) 





459.22 


1.06 


133.64 


1.02 


2 


459.38 


1.11 


133.88 


1.12 


3 


459.43 


1.13 


133.90 


1.11 


4 


459.43 


1.14 


133.90 


1.13 


5 


459.43 


1.13 


133.86 


1.15 


6 


459.45 


1.14 


133.86 


1.16 


20 


459.49 


1.11 


133.83 


1.18 



LAO core levels 

n La 4d 5/2 BE La 4d 5/2 FWHM Al 2s BE Al 2s FWHM 
(eV) (eV) (eV) (cV) 

2 103.20 1.84 119.55 2.00 

3 103.18 1.86 119.53 1.98 

4 103.17 1.83 119.52 1.99 

5 103.04 1.89 119.39 2.08 

6 102.86 1.88 119.24 2.06 
20 103.10 1.81 119.48 1.97 



TABLE I. Ti 2p 3 / 2 , Sr 3ds/2 La 4ds/2 and Al 2s binding energies (BE) and full width at half maximum (FWHM) as a function 
of LAO layer thickness. All binding energy values shown have an error bar of ±50 meV and all widths have an error bar of 
±70 meV. All data shown was acquired from spectra measured using hv — 3 keV and the total energy resolution was 500 meV. 



these electrons before the critical LAO thickness for con- 
duction is reached is still unexplained. 29 Defects, such as 
oxygen vacancies or adsorbates 9 could be the source of 
these electrons, and the occupation of Ti 3d states for 
n < 4 seen here and in Refs. [9 and 29] point to a tran- 
sition from localized to delocalized behavior for at least 
partially pre-existing Ti 3d electrons at < n < 4, rather 
than a transition from zero to finite Ti 3d occupation at 
the critical LAO thickness. Seeing as an x-ray beam is 
an integral part of all core level photoemission exper- 
iments, photo-induced electrons are another possibility 
that need consideration, 29 and we will return to these is- 
sues in Section IV. Naturally, the presence of Ti 3+ ions 
in these systems for n > 4 is in accordance with the elec- 
tron transfer model, a model which also predicts large 
core level shifts in the LAO layers. Thus, in order to 
investigate the presence of these shifts we now move on 
to a discussion of the La 4d and Al 2s core levels for all 
samples, using the Sr 3d level of the STO substrate of 
the same sample as an internal energy reference. 

Fig. 2a shows typical La 4d and Sr 3d spectra for 
n = 2 and n — 6. The Sr 3d spectrum for STO (n = 0) is 
shown for comparison. As expected, the energy splitting 
resulting from spin-orbit interaction (Sr 3d 5 / 2 -3d 3 / 2 , La 
4d5/ 2 -4d 3 / 2 ) is unchanged as the LAO overlayer thickness 
increases. In contrast to the STO core levels for n > 2, 
however, there is a clear shift for the La 4d peak with 
increasing LAO thickness. In the inset to Fig. 2a, the La 



4d5/2 peak maximum is shown for all samples (n = 2 — 6), 
showing the core level shift to lower binding energy in all 
cases, whereby for n > 4 this shift is the most striking. 
The total shift between n = 2 and n — 6 amounts to 

- 300 meV. 

Thus, what Fig. 2a shows is a clear thickness- 
dependent shift of the LAO core levels to lower binding 
energy between n — 2 and n = 6, that on a qualita- 
tive level agrees with the electronic reconstruction model. 
Having said that, it is also immediately evident that the 
La 4d core level binding energy changes by much less 
than the simple expectation of 0.9 eV per LAO unit cell. 
In Fig. 2b we show the La 4d spectrum simulated for 
ri = 3, assuming validity of the simple electronic recon- 
struction model. We use the La 4d HAXPES spectrum 
from La 2 _ 2x Sri +2a; Mn 2 07 with x = 0.425 as a model 
lineshape, as we know there is no internal potential build- 
up in this system 32 . For n = 3, the levels in each of the 
LAO unit cell layers are expected to be shifted by ~ 0.9 
cV with respect to the LAO unit cell layer below. 7 Con- 
sequently, we sum up shifted contributions representing 
emission from each LAO unit cell layer, and then take an 
average weighted by e~ d / A , where d is the depth of the 
layer and A the electron inelastic mean free path of 4 nm 
in this case, as this weighted average better represents 
what could be expected to be measured in experiment. 

It is clear that the simulation in Fig. 2b is too broad 

- visible for example in the almost complete filling in of 
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FIG. 2. (color online) a) La 4d and Sr 3d core level photoe- 
mission spectra for selected STO/(LAO) n samples. The Sr 3d 
data were measured with hv — 2150 eV, the La 4d data with 
hv = 3000 eV. The inset shows a clear shift of the La 4d maxi- 
mum close to 103 eV to lower binding energies with increasing 
LAO layer thickness, whereas the Sr 3d spectra essentially do 
not shift in energy between n = 2 and n = 6. b) A simulated 
n — 3 spectrum for the La 4d core levels which are shifted by 
0.9 eV per LAO unit cell. The colored, dashed spectra are the 
layer-resolved spectra for the first (interface side, red), second 
(orange) and third (surface side, green) LAO unit cell, whose 
A-weighted sum then gives the total (black) spectrum, c) The 
La 4d 5 / 2 peak FWHM as a function of LAO layer thickness. 



the minimum between the two main La 4d peaks at 105 
eV Eb ■ The experimentally observed width is essentially 
unchanged for n > 2, as is shown in Fig. 2c and Table I 33 , 
and there is no indication that the core level peaks consist 
of a sum of more than one atomic core level spectrum. 
The data of Fig. 2 show that although there is an energy 
shift of equal sign to that expected from the electron 
transfer picture, it differs from the expectations of the 
simple model in two important respects: 

• we observe a rigid shift rather than a layer-by-layer 
potential build-up 

• the shift's magnitude is a (small) fraction of the 
expected value of n times 0.9 eV/layer from the 
electronic reconstruction model. 
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FIG. 3. (color online) Evolution of the energetics in the 
LaA103 overlayer tracked using the La 4d and Al 2s core 
level binding energies vs. LaA103 thickness, n. a) Data up 
to n = 6. To indicate the spread in energy shifts encoun- 
tered, in addition to the main data-points (see legend), the 
smaller grey symbols show energy shift data from all experi- 
ments we conducted on STO/LAO interface samples, includ- 
ing using Al:K a radiation, and using samples from other thin 
film growth batches. For comparison with experiment, one 
fourth of the expected core level shift based on DFT calcula- 
tions for a vacancy free interface is also shown as an average 
weighted by the inelastic mean- free path length in HAXPES. 
The results of these calculations are also shown in Table II. 
b) Energy shift data including n = 20. 



Fig. 3 shows the observed core level shifts for all sam- 
ples discussed here, presented in terms of the binding 
energy difference between the La or Al core level con- 
cerned and the Sr core level from the substrate of the 
same sample: A£ B (ii)=E B( „ )Sr3d) - E B { n ,x) with x = 
La 4d or Al 2s. In order to be able to plot both the La 
4d and Al 2s core level shifts (relative to those of the 
substrate) on the same scale, the extrapolated value of 



AEg(n) to n = 0, which we denote AEb has been sub- 
tracted from A_Eg(n). The total shift between n = 2 and 
n = 6 amounts to ~ 300 meV. Also shown in Fig. 3a is an 
e~ d / A -weighted average of the core level shifts from our 
DFT calculations for a vacancy free STO/LAO interface, 
the results of which are shown in Table II. The details 
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of these calculations are given in Section VI. In order to 
be able to compare to the results of earlier theoretical 
predictions, 7 we scaled the theoretical core level shifts 
by a factor of E g /E® FT = 3.2/2 = 1.6 with E g being 
the experimental band gap of STO and E^ FT the DFT 
value of the STO band gap. Subsequently the data were 
divided by a factor of four, to fit on the plot shown in Fig. 
3a. We note here that a simple '0.9 eV energy shift per- 
unit cell' approximation gives a similar result to the data 
shown here. The comparison between theory and exper- 
iment illustrate two things: (i) that the n-dependence of 
the core level shifts observed in experiment is different to 
the DFT prediction, where in the latter stronger shifts 
are expected for n < 4 and essentially no shifts for n > 4; 
(ii) the experimental energy-shift data lie well below one 
fourth of the values predicted by the charge transfer pic- 
ture. Thus, real-life samples lack both large core level 
shifts and, as seen in Fig. 2, the connected core level 
broadening. This forces us to consider that the simple 
layer-by-layer potential build-up picture is not captur- 
ing the main physics at work in real systems. Taken at 
face value, the severe reduction of the potential build-up 
manifested by the reduced core-level shifts indicates that 
other effects, such as defects or reconstructions, should 
be taken into account. 

In Ref. 8, XPS data from MBE grown samples 
recorded using Mg: K a radiation hv = 1253 eV were 
reported. There, a total shift ~ 180 meV was observed 
between samples with n — 2 and n — 6 unit cells of 
LAO. Not only are these shifts 60% of those we observe, 
the sign of the shift is opposite to that clearly visible in 
our raw data. The MBE samples used in Ref. 8 were 
grown in an oxygen partial pressure of 10~ 7 Torr, so it is 
certainly conceivable that they contain significantly more 
oxygen vacancies compared to the 'high pressure' sam- 
ples studied here (po 2 = 2 x 10~ 3 mbar). We will return 
to this point later on. 

The energy offset we observe between n = 4 and n = 6 
of ~ 300 meV (see Fig. 3) is of the same sign and com- 
parable magnitude to that inferred from recent tunneling 
transport data on low O2 partial pressure PLD-grown 
and post-annealed films 34 . These were observed to en- 
ter a Zener tunneling regime at nominally 18.5 unit cells 
of LAO, which would mean ~ 300 meV per unit cell 34 . 
Taken at face value, our core level shift data for n < 6 
could be argued to be consistent with a picture in which 
the critical thickness at which the bottom of the STO 
conduction band is crossed by the top of the LAO va- 
lence band would be much greater than n = 6. A naive 
extrapolation of our data gives a critical thickness as high 
as 35 UC. In order to clarify this we measured an n = 20 
sample, whose core level shift data are shown in Fig. 3b. 
Just as for the n < 6 data, there is no sign of signif- 
icant broadening of the La and Sr core level lines for 
n = 20, ruling out a unit-cell- by-unit-cell energy shift. 
Furthermore, the binding energy shift for n — 20 is in 
fact comparable to that of the n = 5 sample, illustrating 
that the LAO energetics probed via the core levels do not 



simply saturate at some value, but that the behaviour is 
more complex. In any case, as our experimental data 
for n = 20 are inconsistent with the existence of a sec- 
ond critical thickness at high n and with the presence 
of a step-wise unit-cell by unit-cell potential build-up in 
the LaA103, it would seem unlikely that the picture of 
an electronic reconstruction at very high n is relevant in 
these samples. 



IV. PHOTON DOPING 

At this point it is clear that our data do not display 
the results of a major unit-cell by unit-cell energy shift. 
The question arises whether the nature of the experi- 
ment itself influences the electrostatics of the sample. 
The incident hard x-rays used in the experiment create 
electron-hole pairs that could redistribute in the LAO so 
as to cancel the potential build-up and therefore flatten 
the bands. Indeed, it has been observed that any light 
with an energy greater than the STO band gap dramati- 
cally reduces the resistivity of these systems in transport 
experiments. 35-37 

In the case that x-ray-induced carriers do play a role, 
they should not only reduce the potential build-up in 
the LAO overlayer, but should also neutralize the band 
bending or confining potential for the carriers in the STO 
substrate. Therefore, if a sufficient number of electron- 
hole pairs were to be doped into the system in this way, 
both the potential build-up in the LAO and the band 
bending in the STO should disappear entirely. In such a 
scenario - without the confining potential - the Ti 3d elec- 
trons present in the LAO/STO system are then free to 
diffuse through the millimeters thick substrate. Thus, if 
the system is fully dominated by x-ray-induced electron- 
hole pair creation, one would expect both no potential 
build-up in the LAO and no spatial confinement of the 
Ti 3+ ions. 

In a perfect (unreconstructed), defect-free 
LaA10s/SrTi03 sample where no photon-doping is 
present at all, the electronic reconstruction scenario 
would lead to the expected strong potential build-up 
in the LAO and to confinement of the Ti 3+ signal to a 
thin layer of order a few unit cells around the interface. 
Naturally, between this one extreme and the other, in 
which photodoping has flattened all potentials and band 
bending, there is a grey area in which both potential 
build-up in the LAO and Ti 3+ confinement both de- 
crease with increasing photon-induced electron-hole pair 
concentration. One remark is that the band bending 
confining the Ti 3d electrons is suggested to be of 
the order of several hundred millivolts, 31 whereas the 
expected potential build-up in the LAO is an order of 
magnitude greater. Because of this one could expect the 
Ti 3+ confinement to be eliminated first as the possible 
photon-induced electron-hole pair concentration were to 
be increased. Upon a further increase of hypothetical 
photo-doping the LAO bands would then flatten further, 
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FIG. 4. (color online) a) The Ti 2p spectrum for the low 
quality substrate sample compared to the n = 5 sample from 
the high quality series. The sample with the lower quality 
substrate shows a Ti 3+ shoulder which is clearly much bigger 
than that of the high quality sample. The spectrum of the 
low quality substrate was shifted in energy to coincide with 
that of the high quality film, b) The Ti 2p spectrum for the 
low quality substrate, n — 5 sample at different excitation en- 
ergies between 1486 eV and 6 keV. These spectra have been 
normalized on the main Ti 2p 3 / 2 peak. The intensity of the 
Ti 3+ (3d 1 ) shoulder shows a clear kinetic energy dependence, 
in turn related to the inelastic mean-free path length, indi- 
cating confinement of the q-2DEG over 2.5nm, as shown in 
panel (c) in which the Ti 3+ /Ti 4+ intensity ratio is plotted as 
a function of A. 



ultimately resulting in the scenario outlined above. 

As discussed above, in the regular samples we see little 
or no angular dependence of the Ti 3+ shoulder (see Figs, 
lb and c), which - taken at face value - would indicate 
a thickness of the quasi-2DEG of > 8 nm. However, a 
second n = 5 sample was grown under the same condi- 
tions as the high quality series, but on a lower quality 
substrate, which exhibited point defects in an AFM in- 
spection (taken after the usual surface termination pro- 
cedure, but prior to the LAO deposition). As Fig. 4a 
shows, the lower quality substrate gives an LAO/STO 
sample with a large Ti 3+ contribution: the shoulder to 
the Ti 2p main line being at least a factor two higher 
than the one in the regular n = 5 sample. By varying 
the kinetic energy of the photoclectrons via changing the 
photon energy step-wise from 1486 eV to 6 keV a depth 
profile of the Ti 3+ /Ti 4+ concentration ratio of the lower 
quality substrate sample was made. The Ti 3+ /Ti 4+ in- 
tensity ratio extracted from panel (b) of Fig. 4 is plotted 
as a function of A in panel (c). For the special case of 
normal emission (6 — 0), we can fit Eq. 1 to the data 
and extract a thickness for the q-2DEG of 2.5 ±1.5 nm, 
with the fit parameters for the best fit and the range of 
parameters that give reasonable fits indicated in the fig- 
ure legend and as an orange band in Fig. 4c. Both the 
fit parameters and their error bars are in the same ball 
park as the ones reported in Ref. 29. 

If catastrophic photon-induced potential flattening in 
the LAO were operative under the conditions used for 
these HAXPES experiments, then the same would pre- 
sumably hold for the confining potential in the STO, and 
thus no confinement of the Ti 3+ to a depth within ca. 
2.5 nm of the interface would be observed. Thus, the fact 
that clear confinement is observed in this particular n = 5 
sample argues against dominance of photon-induced po- 
tential flattening in these experiments. Furthermore, it 
is germane to note that the La 4d and Al 2s core level 
binding energies of the lower quality substrate n = 5 
sample are very similar to those of the n = 5 sample 
from the regular, high quality series (see Fig. 3a). The 
fact that these two n = 5 samples, which were measured 
in the same photon beam, have such similar La 4d and 
Al 2s core level binding energies, yet such different Ti 3+ 
confinement suggests - but does not prove - that photon 
doping is not the cause of the lack of band bending in 
the high quality sample. Which in turn suggests that 
the level of photon-induced band flattening in the high 
quality sample series is not overly large. 

In an experiment in which we held an n — 3 sample in 
the dark for 40 hours and then measured the Sr 3d, Al 
2s and La 4d core levels using A1:K Q radiation at 10% of 
the regular photon flux we observed no change in both 
the core level binding energies and the core level widths 
as the measurement progressed over the timescale from 
30 minutes to several hours in the beam. Also no change 
was observed after 24 hours, which would be twice the 
amount of time required to accumulate 0.5 electrons per 
UC at the STO /LAO interface at this intensity, as we will 



discuss in the following, in which we make 'back of the 
envelope' estimations of an upper bound for the num- 
ber of electron-hole pairs the A1:K Q x-ray source used 
in the lab experiments could generate in the sample. A 
simple derivation shows this upper bound to be equal to 
Eq. 2. Assuming that the primary photoelectron current 
itself corresponds to a negligible amount of the energy 
deposited in the sample in the form of 1486 eV photons, 
we reason that each photon represents enough energy to 
excite across the STO energy gap 1486/3.2 = 464 times 
or 265 times across the LAO energy gap. Given the total 
photon flux (j> of the A1:K Q source of 10 11 photons/s 38 , 
impinging on a sample area A of 5 x 5 mm 2 , and taking 
the x-ray attenuation length in LAO Xlao to be 0.32 



/jm and similarly Xsto = 1-68 ptm 39 we arrive at an area 
density of electron- hole pairs heh in an LAO layer of 
thickness d LAO = 3 UC of 7.4 x 10 10 cm^s -1 . Natu- 
rally, including processes which certainly take place such 
as energy-loss to phonons would reduce the estimate of 
the number of possible gap excitations from this upper 
limit. The number of electron-hole pairs required to can- 
cel the potential build-up is 3.3 x 10 14 cm -2 (0.5 electron 
per UC). Thus, if only electron hole pairs in the LAO 
contribute to cancellation of the potential build-up, the 
system would reach this level after some 75 minutes in 
the ALKq, x-ray beam. Reducing the photon flux by a 
factor of ten (as we did in the lab experiments) would 
then increase this time to a little over twelve hours. 



lao = f ^ X LAO (e- d ^^/ x ^Q - 1) 

lEH A ■ E% AO \ LAO {e- d ^o/\LAo _ i) + X STO (e-dLAo/\sTo) 



(2) 



Now we consider a scenario in which electron hole pairs 
created in the top 16 nm (twice the value we found for 
confinement of the q-2DEG in Fig. lc) of the STO also 
contribute to the band flattening in the LAO layer. The 
additional number of electron-hole pairs created is given 
by Eq. 3, with d* the thickness of the STO layer that 
donates electron-hole pairs. This number amounts to 
1.6 x 10 12 cm _2 s -1 , which means one could reach 0.5 
electrons per UC already after 3 minutes at full intensity 
and 30 minutes at 10% intensity. Given the low photoion- 



isation cross-sections and resulting modest signal to noise 
ratio in these experiments, it is not possible to determine 
whether there are any changes in core level binding ener- 
gies or widths on this time scale. Thus, without depend- 
able information regarding to what extent electron-hole 
pairs in the STO substrate play a role in possible flatten- 
ing of the bands, and how much of an overestimate our 
upper bound represents, these estimations indicate that 
via this route alone we cannot make a clear call as to the 
effect of photon doping in these experiments. 



,STO 
l EH 



<j) ■ hv X S To{e~ d * /XsTO - e - dLAO / XsTO ) 

A ■ E^ TO XsTo{e~ d * /XsTO - e- dLA °/ XsTO ) + X L Ao(e~ dLAo/XLAO - 1) + AsTo(e~ d * /AsTo) 



(3) 



Before concluding this section dealing with aspects of 
the experiments themselves that could affect the observed 
energetics, we briefly discuss sample contacting. 40 As re- 
lated in the experimental section, the surface of the films 
is additionally contacted by means of colloidal graphite, 
leaving a strip of 2 mm free on which the x-ray beam 
impinges. We found this to be an effective manner to 
prevent charging. A recent paper 41 uses DFT techniques 
to study the effect of adding metallic contacts (Ti, Al, 
Pt) on top of LAO on STO. The conclusions are that 
such metallic contacts provide a significant charge trans- 
fer to the Ti states of the STO, turning them metallic and 
neutralizing the electric field in the LAO. For Ti and Al 
these effects are predicted to be largest, as strong cova- 
lent bonds between the outermost AIO2 layer of the LAO 
and the metal result. For Pt, the bonding is weaker and 
the calculations predict a residual, but reduced electric 
field in the LAO. 41 



In our case, a mm away from the region of the films 
measured using HAXPES there is an overlayer of col- 
loidal graphite. This carbon film can be simply removed 
by wiping with a tissue impregnated with isopropanol, 
thus everyday lab practice confirms one's chemical intu- 
ition, expecting no covalent bonding and charge transfer 
between the graphite colloids and the LAO surface. Con- 
sequently, one would not expect this adjacent graphite 
contact layer to affect the energetics observed from the 
exposed LAO surface from which the data presented here 
are reported. Indeed: when contacting a sample using 
a minimal amount of silver conductive paint only on 
the bottom (substrate side) of the sample we observe 
no change in any core level binding energies (except due 
to electrical charging) or peak shape when compared to 
truly contacting the STO/LAO interface of the sample 
using the colloidal graphite. 

Before wrapping up as regards the thickness depen- 
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dent energy shifts observed in LAO on STO from our 
HAXPES data, it is only fair to mention that a black or 
white proof either way on the issue of the possible effect 
of photo-doping or surface photo-voltage effects (and a 
proof that the sample contacting method plays no role 
whatsoever) is not possible at this stage, but our con- 
sidered opinion is that the core-level shifts we report in 
Figs. 2 and 3 and their deviation from the simple 0.9 eV 
unit-cell by unit-cell potential build-up scenario are to be 
taken seriously. 



V. EXPERIMENTAL DATA SUMMARY 

Summarizing the experimental results: 

• We have detected a clear low binding energy shoul- 
der as part of the Ti 2p 3 / 2 core level photoemission 
feature, consistent with Ti 3d states being popu- 
lated with electrons for samples with LaAlC>3 thick- 
nesses n > 2. 

• The relative core level binding energies of the LAO 
layers decrease as a function of LAO thickness in a 
manner consistent with an increase in the LAO po- 
tential. This apparent potential build-up increases 
at a higher rate for n > 4, but decreases again for 
n > 6. 

• Although the sign of the binding energy shifts in 
LAO agree with the simple electrostatic model for 
potential build-up, their magnitude - some ~300 
meV between n = 2 and n = 6 - is much smaller 
than the predicted 0.9 eV per unit cell. 

• In addition, the data do not support the existence 
of a unit-cell by unit-cell alteration of the electro- 
static environment in the LAO overlayer, as the 
LAO core level lines show no measurable increase 
in width as the overlayer thickness increases. 
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FIG. 5. (color online), a) The n-dependence of the potential 
in the LaAlOs expressed via the relative changes in the La 
5s shallow core level binding energy extracted from the DFT 
calculations for vacancy free and oxygen deficient LAO/STO 
systems described in the text. Points describing a positive 
slope reflect a trend to decreasing La core level binding energy, 
b) The 2x2 lateral supercell used in the DFT calculations. 
Oxygen vacancies (25%) were located in the outermost AIO2 
layers, c) Layer-by-layer data from the DFT calculations for 
the system shown in panel (b). The right hand panels show 
the n-dependence of the energetics in the STO substrate and 
the left hand panels the evolution of the energies across the 
LaAlOs films. The yellow shading indicates regions with LAO 
potentials such that the La core levels move to higher bind- 
ing energy, opposite to the direction of the in-built potential 
scenario. 



VI. DFT CALCULATIONS 

As touched upon earlier, oxygen vacancies could form 
part of the explanation of both the Ti 3d occupation for 
n = 2, and the reduced potential build-up compared to 
the simple electron-transfer model. Therefore, we turn 
to an investigation of the effect of oxygen vacancies on 
core level shifts by means of DFT calculations. In the 
computer, we have introduced vacancies on 25% of the 
oxygen sites in the outermost AIO2 layers of the LAO. As 
this is the energetically most favorable position for oxy- 
gen vacancies in the LAO overlayer 12 we will only con- 
sider this location for oxygen vacancies. For discussion of 
oxygen vacancies created in other layers of the system as 
well as in the outermost AIO2 layer, we refer the reader 
to Ref. 11. A 25% oxygen vacancy concentration would 
nominally provide enough carriers (half an electron per 



unit cell) to completely compensate the potential build- 
up caused by the polar catastrophe, eliminating the need 
for electronic reconstruction. In this sense, the calcula- 
tions presented are only a limiting case scenario, but are 
instructive nonetheless. 

For the DFT slab calculations we focus on symmetric 
(ra,ll/2,ra) LAO/STO/LAO slabs in which n unit cells 
of LAO (2n atomic layers) on either side of 11/2 unit 
cells of STO - resulting in two n-type interfaces - are 
repeated periodically, separated by a vacuum region of 
approximately 15 A, see Fig. 5b. Because the samples 
are grown on STO substrates, we fix the in-plane lat- 
tice constant at the calculated equilibrium value of STO. 
We introduce oxygen vacancies on both LAO surfaces 
in such a way that inversion symmetry in the central 
layer is maintained. Oxygen vacancies at the surface of 
the LAO are modelled in a 2x2 lateral supercell and all 
atoms are allowed to relax internally. The initial atomic 
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structure is undistorted, without considering TiOe and 
A1C>6 octahedron tilting and rotation. Due to symmetry 
breaking induced by the oxygen vacancies at the surface 
AIO2 layer, atomic relaxation will lead to an out-of-plane 
tilting of A1C>6 octahedra. 11 Even so, it still cannot de- 
scribe a possible in-plane rotation of GdFe03-like distor- 
tions at the interface, 42 the investigation of which goes 
beyond the scope of this article. The computational de- 
tails have been given elsewhere. 12,42 The effect of strong 
correlations, which have been shown to be important in 
STO/RO (R = La, Pr, Nd, Sm) interfaces, 43 are not 
treated explicitly here. Although their inclusion in terms 
of LDA+U may modify the electrostatic potential pro- 
file slightly, it would not change our main theoretical 
conclusions, which are dominated by the polarity of the 
system. 12 The La 4d and Sr 3d core levels convenient for 
experiment are too deep in binding energy for the cal- 
culations, so instead we use the shallower La 5s and Sr 
4s core states to track the potential changes at the La 
and Sr sites. In Figs. 5b and 5c we introduce a unit 
cell index, m, which is for STO at the LAO interface, 
negative going deeper into the substrate and m > 1 in 
the LAO. We calculated the density of states (DOS) for 
samples with LAO layer thicknesses 1 < n < 6, projected 
onto the last few layers of the STO substrate and each of 
the LAO layers on top. From these data the La 5s and 
Sr 4s core level binding energies could be extracted, and 
these are shown in Table II for the 25% oxygen vacancies 
case and for the vacancy free case. The internal energy 
reference used for both core level energies is the Sr 4s 
core level in the central SrO layer of the STO substrate 
furthest from the interface with LAO (to = —2). 
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TABLE II. The binding energy difference between Sr 4s (m = 
—2) and La 5s in eV as calculated for both the 25% oxygen 
vacancy case and the vacancy free case as a function of LAO 
thickness in UC for all values of m. 



The upper (red) curve in Fig. 5a shows the theory data 
for the mean-free-path weighted average of the change in 
the binding energy of the La 5s core level (with respect to 
an extrapolated value for 77 = 0) as a function of the LAO 
layer thickness (also shown in Fig. 3a, but without the 
3.2/2 bandgap correction). The curve shows a A weighted 
energy shift equivalent to 0.3 eV per 77 up to 77 = 4 and 
is then roughly constant thereafter. The per-layer shifts 
shown in Table II for the vacancy free case correspond, 
for instance for the 71 = 5 sample, to an average potential 
build-up of about 0.5 eV per layer and are in agreement 
with the results of Ref. 44. DFT calculations underesti- 
mate the magnitude of the STO band gap, putting it at 
2.0 eV, whereas the experimental value is 3.2 eV. 44 This 
is why - within the electronic reconstruction model - the 
core levels are in reality expected to shift by 0.9 eV per 
unit cell, rather than the 0.5 eV per unit cell the DFT 
data indicate. 

Before discussing the lower (green) curve in Fig. 5a, 
we turn to the layer-resolved core level data from the 
calculations in Fig. 5c. Such layer-resolved data are in- 
accessible in a HAXPES experiment, but nonetheless the 
discussion of this theory data can give some insight into 
the more detailed processes underlying the type of global, 
77-dependent energy shifts detectable in experiment. 

The right hand panels of Fig. 5c indicate that the 
(near)-interfacial Sr 4s core level binding energy is greater 
than that in the 'bulk' of the STO substrate (to = —2) 
by about 200 meV. This shift is a consequence of the 
existence of an LaO + layer immediately atop the Ti02- 
terminated STO, and the theory data show it to be prac- 
tically independent of LAO thickness 77. This, in fact, 
agrees with the experimental situation for the Sr 3d and 
Ti 2p core levels, as described at the outset of Section III 
and shown in Table I. 

The left hand panels of Fig. 5c illustrate the devel- 
opment of the La 5s binding energies (referenced to the 
Sr 4s of the substrate [m = —2]) for the system with 
25% O vacancies. Looking first at the interface (to = 1), 
one sees that the binding energy of the La 5s states in- 
creases steadily as the LAO thickness increases, leading 
to a reduction of BEs^-BE^^s. On the other hand, 
the STO-referenced La 5s binding energy at the LAO 
surface (i.e. for m — n) remains relatively unchanged as 
n increases. 

The A-weighted sum of the unit-cell by unit-cell 
changes shown in Fig. 5c gives an overall trend shown in 
the lower (green) data-points in Fig. 5a. Here, in analogy 
to the data plotted in Fig. 3 we consider the 77-dependent 
change in La shallow core level binding energy compared 
to the hypothetical value for n — 0, the latter derived 
from an extrapolation of the n = 6 — > 1 values through 
to 77 = (a quantity referred to as A-Eb in the ordinate 
label). Fig. 5a shows that, for the case of 25% O vacan- 
cies, the La 5s core levels (referenced to a 'bulk' Sr 4s 
level in the center of our slab) shift towards more nega- 
tive values, indicating a relative increase in the La 5s core 
level binding energy as 77 increases for the case of 25% 
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O vacancies 45 . This shift is opposite to the vacancy free 
case (no oxygen vacancies, red symbols in Fig. 5a), and 
illustrates how the core level shift to lower binding ener- 
gies expected from the electronic reconstruction scenario 
can be negated and even turned around by the presence 
of oxygen vacancies in the LaAlC>3. Because of the non- 
trivial layer-by-layer changes shown in Fig. 5c, LAO core 
level binding energies for oxygen vacancy concentrations 
between zero and 25% cannot simply be arrived at by ap- 
propriate interpolation between the red and green points 
in Fig. 3a. Nevertheless, the 'limiting-case' presented 
of 25% O vacancies does illustrate effectively that such 
defects - and their electronic consequences - will have a 
strong effect on the energetics in the LAO layers, and are 
able to neutralize or even reverse the potential build up 
in the LaA103 layer. 

This effect could possibly explain the results of Ref. 8 
as regards the opposite n-dependence of the LAO core 
level binding energies compared to the data presented 
here, as the samples in Ref. 8 were grown in a low 
O2 partial pressure. The experimental data would be 
consistent with our moderately high O2 partial pressure 
PLD-grown films being in between the two extremes of 
the defect-free, pure charge transfer picture and a high 
oxygen defect scenario. Interpolating between the results 
from the defect free and 25% oxygen vacancy calculations 
may suggest that the observed energy shifts are compat- 
ible with an oxygen defect concentration in the top AIO2 
layer of our samples of 12-16%. However, we emphasize 
that the experimental data show no sign of a unit-cell- 
by-unit-cell energy offset, whereas - for less than 25% O 
defects at the LAO surface - the DFT calculations would 
predict such a step- wise shift. Because oxygen vacancy 
formation is driven by the incipient polar catastrophe, 
this driving force will be diminished upon the formation 
of vacancies, and the 25% concentration considered here 
will be a limiting case. In reality, it is therefore not un- 
reasonable to find a lower oxygen defect concentration 
and -in addition- other kinds of (charged) defects or non- 
stoichiometries that also play a role in the flattening out 
of the built-in potential shifts in real samples. 

VII. CONCLUSION 

In conclusion, we have studied STO /LAO heterointer- 
faces, grown under moderately high oxygen partial pres- 
sure PLD using hard x-ray photoelectron spectroscopy. 
The aim is the investigation of possible built-in fields in 
the LaA103 by tracking the binding energy differences 
between the STO and LAO core levels as a function of 
LAO thicknesses from n — 2 to 6 unit cells. 

Our data show a positive core level shift with increas- 
ing LAO thickness (i.e. binding energies of the LAO core 
levels decrease) : the sign of this shift meets the expecta- 
tions from simple built-in potential considerations. The 
magnitude of the binding energy shift, however, is very 



low: only ~ 300 meV, the majority of which occurs be- 
tween n = 4 and n = 6, and is not of the unit cell by 
unit cell type. The observed shifts are much lower than 
the 900 meV per unit cell expected in the simplest ionic 
polarization models, and remain so even for the large n 
value of 20. 

Our DFT calculations show that oxygen defects in the 
AIO2 surface layer of the LAO offer one possible expla- 
nation for this reduction of the internal potential, as 25% 
oxygen defects (sufficient, in principle, to deliver 0.5 elec- 
trons per Ti ion at the interfacial layer) are sufficient to 
even reverse the direction of the energy shift of the La 
levels in there LAO. 

The facts that we experimentally observe occupation 
of the Ti 3d states as a shoulder to the Ti 2p signal al- 
ready for n = 2 and n = 3, and that there are only very 
modest energy shifts observed per additional LAO unit 
cell strongly suggest that the critical LAO thickness of 
n = A seen in transport studies is not directly related to a 
built-in potential-induced shift of the LAO valence band 
to meet the STO conduction band states. Consequently, 
it is apparent that in real materials, the electronic recon- 
struction scenario is pre-empted by other effects, of which 
we show oxygen vacancies to be one plausible possibility. 
Influence from photo-induced charge carriers cannot be 
excluded in our photoionization-based experiments, but 
we consider it unlikely that these effects dominate the 
physics we observe. 

Finally, the spectroscopic data do not single out n = 
4 as being special, in contrast to transport data from 
numerous different groups. This suggests that although 
the unit-cell by unit-cell potential shift does not - in the 
end - appear to be the scenario met in real materials, 
there still must be a mechanism providing 'action at a 
distance', for example from the outermost electrostatic 
termination layer of the STO/LAO film stack, feeding 
back to the localized or itinerant character of the states 
at or near to the STO interface, in this way generating 
the observed critical thickness seen in transport. 
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